End of multi-field infiation and the perturbation spectrum 



Jinn-Ouk Gong* 



International Center for Astrophysics 
Korea Astronomy and Space Science Institute 
Daejeon, Republic of Korea 

February 2, 2008 



Abstract 

We investigate the dynamics of inflation models driven by multiple, decoupled scalar 
fields and calculate the Hubble parameter and the amplitude of the lightest field at the end 
of inflation which may be responsible for interesting, or possibly dangerous cosmological 
consequences after inflation. The results are very simple and similar to those of the single 
field inflation, mainly depending on the underlying spectrum of the masses. The mass 
distribution is heavily constrained by the power spectrum of density perturbations V and 
the spectral index ng. The overall mass scale gives the amplitude of V, and ng is affected 
by the number of fields and the spacing between masses in the distribution. The drop-out 
effect of the massive fields makes the perturbation spectrum typically redder than the 
single field inflation spectrum. We illustrate this using two different mass distributions. 
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1 Introduction 



After its advent, inflation [Ij has been very successful to solve many cosmo logical problems of the 
standard hot big bang cosmology and its predictions are consistent with the recent observations 
[21 E]: inflation can explain the homogeneity, isotropy and flatness of the observable universe, 
as well as the origin of the large scale structure in the universe. Interestingly, the simplest 
possibility, i.e. the model of chaotic inflation with quadratic potential V{(f)) = m^0^/2 is still 
consistent with the WMAP three year results [3] . Thus it is natural to try to build such a fully 
realistic inflationary scenario in the context of string theory, the most promising candidate of 
the fundamental theory of nature. However, it is still unclear how to flnd successful inflationary 
scenarios in string theory. Especially, if the vacuum expectation value of the inflaton fleld is 
beyond the Planck scale mpi, which is the case for the simplest m^0^ chaotic inflation, generally 
large radiative corrections are expected and they would modify the inflaton potential V{(j)), and 
the predictions based on the original form of the effective potential will be no more reliable. 

Fortunately, after the de Sitter vacua construction in type IIB theory [1], there have been 
considerable advances to this end for recent years [5]. An interesting scenario named N-flation 
was recently proposed [6] where the m^0^ inflation model is achieved by the assisted inflation 
mechanism [7]: many (A^ 3> 1) flelds combine to drive inflation with the displacements being 
less than mpi. These flelds are considered to be the string axion flelds, which are known 
to exist abundantly in string theory, and they have flat potentials even after all the moduli 
are stabilised. In Ref. [S] the special case where all the masses are the same is highlighted, 
and subsequently more general extensions were investigated in Refs. [H [9] . The general and 
qualitative evolution of the flelds is however clear: the most massive flelds roll down their 
potentials to the minima, and the light flelds are responsible for the later inflationary stages. 
Here one important thing is that, at the end of inflation, not all the flelds are relaxed to the 
minima. For the single fleld inflation with V{(f)) = m^0^/2, as we will see in the next section, 
inflation ends well before reaches the minimum of the potential. Likewise, inflation proceeds 
no more although there are still some flelds whose amplitudes are non-zero. These non-zero 
flelds are of cosmological importance after inflation: they may lead to interesting, or possibly 
dangerous cosmological consequences. For example, they may oscillate and decay into dark 
matter particles [10], they can produce additional perturbations via curvaton mechanism^ |12j . 
and they would be responsible for the preheating [13] of the universe after inflation. They can, 
of course, produce gravitino at an unacceptable level [H], spoiling the successful predictions 
of the big bang nucleosynthesis [15]. The dynamics of the flelds depends on the underlying 
spectrum of masses, which is revealed via the power spectrum V and the spectral index Ug'- 
they depend on the detail of the mass distribution. 

Our purpose in this paper is to study the dynamics of this class of multiple fleld inflation to 
calculate several interesting physical quantities at the end of inflationary regime in connection 
with the power spectrum V and the spectral index rig. Since the lightest fleld keeps the largest 
amplitude at the end of inflation, we will mainly focus on this fleld and its mass: also we can 
describe the evolution of the lightest fleld under slow-roll approximation, and this helps to 
simplify the calculations. This paper is outlined as follows. In Section [21 we flrst briefly review 
the simplest chaotic inflation model with V{(j)) = m'^cfp' /2 and calculate some quantities of the 

"'^Note that for the curvaton mechanism to work properly, we must put more constraints on the parameter 
space of the model under consideration [llj in addition to the conditions for successful inflationary phase. 
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model. In Section [31 we calculate the Hubble parameter, the number of fields with non-zero 
amplitudes and the amplitude of the lightest field at the end of infiation. In Section H] we apply 
the results of the previous section to two interesting mass distributions, and compare with the 
single field infiation case. In Section [5] we summarise and conclude. 



2 Single field infiation 

In this section, we recall the inflation model where accelerated expansion is driven by a single 
inflaton field with simple quadratic potential [1]. The infiaton potential is given by 

l^(</)) = ^mV, (1) 

and with the slow-roll equation of motion for 0, 

3H(j) + m^<\) = , (2) 

the total number of e-folds is ^ 

M= I Edt^ . (3) 

Then, to obtain M > 60, we need 0j > 4\/T5mpi ~ 15mpi: the initial value of the infiaton is 
required to be much larger than m^\. Now, from the acceleration equation 



a 6mpj 

(4) 

3mpj 

we find that infiation ends when 

02>^mV, (5) 

i.e. the kinetic energy of the field becomes greater than half of the potential energy. Using 
Eq. ([2]) to eliminate 0, the value of the Hubble parameter at the end of infiation is estimated 
as 

^^end = -^"^ • (6) 

Note that Eq. ^ is equivalent to the condition that the slow-roll parameter 

becomes larger than 1. Substituting Eq. (EI) into Eq. (JTj) and using Eq. ([5]), we find that at the 
end of infiation has the amplitude 

2 

0end = "^Pl ■ (8) 
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The mass of the inflaton field is heavily constrained by the amplitude of the primordial density 
perturbation^. The leading-ordei^l power spectrum of primordial density perturbations is given 

by [1] 

(9) 



V 

and the corresponding spectral index by 

n.s - 1 



967r^mp[ ' 



5m 



PI 



(10) 



Using the COBE observation 5h = 2VP/5 ~ 1.91 x 10-^ we find 

m ~ 6.07 X lO^^mpi . 
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Note that ris is independent of m and is given by Us — 0.967 when = 4v 15mpi, i.e. 60 e-folds 
before the end of infiatioE0. In Fig. [H we illustrate some results of numerical calculations, and 
in Table [1] we compare the corresponding analytic and numerical estimates. 





Af 


0cnd 




V 




single field 
infiation 


analytic 


100.00 


0.942809 


0.471405 X 10-^ 


6.18101 X 10-^ 


0.966942 


numerical 


100.97 


1.00934 


0.574670 X 10-^ 


5.96853 X 10-'^ 


0.966359 


exponential 
distribution 


analytic 


100.00 


0.121716 


0.101500 X 10-^ 


2.47245 X 10-9 


0.936631 


numerical 


101.43 


0.225956 


0.124816 X 10^^ 


Mar cenko- Past ur 
distribution 


analytic 


99.99 


0.175075 


0.141009 X 10"^ 


4.40243 X 10'^ 


0.955755 


numerical 


101.27 


0.324082 


0.190495 X 10-^ 



Table 1: Analytic and numerical estimates of several quantities of the single and multiple field 
infiation models under the mass distributions, Eqs. (130|) and (!32l) . The quantities for the single 
field infiation are from Fig. [H Here (pend and Hend are written in the unit of mpi, and V and Hs 
are calculated at 60 e-folds before the end of infiation. Note that since it is not easy to derive 
analytic estimates for V and Us for multi-field cases, we only show their numerical results. 
Generally the analytic results agree with the corresponding numerical ones. 



3 Multiple field case 

Now let us consider the infiation case with multiple, uncoupled fields [6]. Here infiation is 
driven by a number of non-interacting string axion fields. If each axion field is displaced not 

^Note that the inflationary energy scale is constrained from the detailed calculations on the amplitude of the 
power spectrum of the primordial gravitational waves [16j . 

■^For the calculations with higher-order corrections in the slow-roll parameters, see e.g. |17) . 

'^Note that, in fact, there exist some level of uncertainty on exactly when the observable large scale pertur- 
bations were generated. While for most models of inflation the corresponding e-fold lies between 50 and 60 [18], 
the values of V and especially Ug are similar: we have Us = 0.963280 and 0.959581 for 55 and 50 e-folds before 
the end of inflation, respectively. Thus, in this paper we keep evaluating V and rig before 60 e-folds. 
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Figure 1: Numerical results of the single field inflation with the initial amplitude 0j = 20mpi 
and the mass m = 10~^mpi. The dotted lines indicate the end of inflation: in the plot of H 
(upper right) the horizontal dotted line is -ffend, given by Eq. ([6]). In the lower row, we set 
A/" = as the end of inflation. 



too far from the minimum, we can write the potential as 

N 

2 



^ 1 

V = Y.\<^h (12) 



1=1 



where we assume that there are N axion fields. The total number of e-folds is then given by 

^= ^yi , (13) 

where (pf'^ is the initial value of (pi. From Eq. (jlj), we obtain 

N 

a 



a 3m^, . 

1=1 
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thus inflation ends when 

N N 



i=l i=l 

Now, using the slow-roll equation of motion for 0j, 

3if0i + m^0i = , (16) 
we can find that under the slow-roll approximation different fields (pi and (pj satisfy the relation 

m 

(pi mf(pi 



thus we obtain 



(17) 



Therefore, we can write (pi and (pi in terms of the lightest field (pN and (pN 



9 / 9 



/,(o) 



/ (p^ 



my V 0(0) / \ JO) 



(19) 



respectively. Substituting Eq. (fT9|) into Eq. (ITSl) and using the slow-roll equation of motion for 
0iv, we find that inflation ends when 



Because (p^^^ /(p^^ < 1 at the end of inflation and rrii > itln always, the last term of each sum 
gives the most significant contribution for that sum. Therefore, the Hubble parameter at the 
end of inflation is 

where ttin is the mass of the lightest field. Interestingly, this is exactly the same form as the 
single field inflation case, Eq. ([6]). Note that at this point, combining Eq. ( |2T1) and the slow-roll 
equation of 0Ar, the kinetic energy of the field 07v is half of the potential energy, 

Kn = ^Vn , (22) 
as the case of the single field inflation. Likewise for (pi we can derive a similar relation 



Now, assuming that total n light fields (n < A^) contribute at the end of infiation, we write 



H. 



2 2 
end — 



[{Kn-u+i + Vn-u+i) + --- + iKN + Vn)] 



(24) 



As can be seen from Eq. ( |23|) . as the mass becomes heavier, the kinetic energy is increasing 
while the potential energy is decreasing. Thus, we can assume that the energy contribution of 
each field roughly coincides. Then, Eq. becomes 



n 



3mpj 



(Kr, + Vj 



N 



(25) 



hence using Eq. f l22|) the amplitude of the lightest field at the end of inflation is given by 



iend 



mpi 



n 



(26) 



which is again similar to the single field case, Eq. ([8]), divided by some numerical factor n. Note 
that this number would be proportional to the total number of fields N unless the masses are 
anomalously distributee^, hence 0^^^ decreases as increases: with larger number of fields, 
the infiationary phase lasts longer until 0Ar reaches smaller value. To estimate the mass of the 
most massive field which contributes at the end of infiation, first we observe that heavy fields 
drop out at the early stages of infiation while light fields drive infiation at late stages. Hence, it 
is reasonable to assume that the fields which are responsible for the infiation of the last e-folds 
have masses within far less than an order of magnitude, otherwise some of them would have 
already quitted the infiationary regime. Thus, we can set the mass of the heaviest field to be 



mN-n+l 



V2: 



TUN . 



(27) 



This number n is dependent on the specific mass distribution as we will see in the next sec- 
tion. We cannot, however, follow the same steps to calculate the power spectrum of curvature 
perturbation^ and the spectral index: following 5J\f formalism [22], they are written as ^ ^ 



V 



E^Y.2 



(28) 



and 



rio — 1 = —4m 



PI 



4 



+ 



(29) 



^If the axion decay constant could be made far larger than mpi, we can obtain consistent infiationary scenario 
with very Umited number of fields, e.g. only two fields [20\. 

^Generally, it is expected that orthogonal isocurvature perturbations are also generated. See, e.g. [21] for a 
review. 
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respectively. From Eq. fl2^ . we can see that the spectrum is always redder than the spectrum 
of the single field infiatioiu. At the end of infiation, the amplitude oi (j)N has significantly 
decreased from its initial value 0^^"*, and thus we could ignore all the contributions except the 
one by (pN in Eq- (|20ll . But before 60 e-folds where the COBE observation is made, 0Ar is 
almost frozen and hence the contributions other than 07v should be taken into account. 

4 Mass distributions 

In this section we illustrate the results of the previous section using two different mass distribu- 
tions of the axion fields. Here we take (pf^ = ■ ■ ■ = = rnp\ for numerical calculations: when 
building infiationary models with masses being given by the underlying physics, the one with 
the largest possible number of e-folds is favoured because it occupies the greatest volume in the 
context of eternally infiating universe [2lj and is preferred a posteriori. This is implemented 
by assigning the maximum initial displacement to the fields with a given spectrum of masses. 

4.1 Exponential distribution 

We first consider the exponential mass spectrum where the fields are distributed uniform on 
logarithmic scale [6l [9] , 

= mje-^'-'^^^ , (30) 

where i = 1, ■ ■ ■ , is the index of the fields, and a is the density of fields per logarithmic 
interval, a is determined once we decide the heaviest mass mi, the lightest mass m^r and the 
number of fields A^. 
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Figure 2: The evolution of the fields. In the left panel, masses are distributed according to 
Eq. ( !30l) where the most massive mass mi = mpi, while in the right panel Eq. (|32l) is applied 
with the average mass m = mpi and the parameter (3=1/2. As in Fig. [H the vertical dotted 
line denotes the end of infiation. 

^Note that this is true not only for the quadratic potential given by Eq. (fT2ll , but for a set of general power 
law potentials [25] , 
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We have performed numerical calculations under this mass distribution with the lightest 
mass being 1/100 of the heaviest mass, i.e. ttln = 10~^mi. We have tried = 250, 500 and 
1000, with mi = mpi, 10~^mpi, 10~^mpi and 10~^mpi for each case. The evolution of several 
light fields is shown in the left panel of Fig. [2l Here, we show only the case where we have set 
mi = mpi and = 250, since for other cases the fields follow exactly the same evolutions except 
for time scale: e.g. it takes 10^ times longer when mi becomes 10^^ times lighter. The dotted 
line indicates the end of inflation. As can be seen, heavy fields drop out of the inflationary 
regime at early times and only a limited number of light fields keep non-zero amplitudes at the 
end of inflation. To estimate the number of the fields which contribute at the end of inflation, 
we can combine Eqs. flTTI) and (1301) to obtain 

(A^-l)log2 , , 

4 log 10 - 

which gives n = 19, 38 and 76 for A^ = 250, 500 and 1000, respectively. Also note that because 
of the relatively large H, many fields undergo dissipation [B] at the early stage of inflation, and 
do not oscillate near the minima. 

In the left panel of Fig. [3] we show the power spectrum V given by Eq. ( l28l) for this case. 
Now, we obtain different spectra for different cases: the mass gives overall amplitude, and the 
number of fields controls the shape of the spectrum, as we will see soon. As can be read from 
Eq. (!28|) . the amplitude of the spectrum is proportional to the sum of the mass squared of 
every field. We can also understand the reason why the spectrum becomes flatter with larger 
number of fields as follows: since here we take (pf^ = mpi for all i, massive fields dominate the 
total energy density at early times. Therefore as these fields are completely dropping out of 
the inflationary regime, the Hubble parameter H is decreasing steeply, making the amplitude 
V decline quickly. This "drop-out" effect of heavy fields weakens when there are many other 
light fields so the contributions by these fields keep the energy density and hence H smoothly 
varying within the slow-roll regime. This is independent of the overall mass scale but dependent 
only on the number of fields and the spacing between masses. Thus, all in all, we can observe 
the same shape of the spectrum with the same number of fields. 

Now we turn to the spectral index Ug. As we have discussed above, here Us is independent 
of the overall mass scale, so in the left panel of Fig. H] we show different indices on different 
A^ only. As anticipated from the discussion above, we obtain larger tilt with smaller number 
of fields. One interesting thing is the initial steep decrease of the index for A^ = 250. This is 
due to the drop-out effect we discussed above: massive fields which occupy a large fraction of 
the energy density drop out of the inflationary regime, so H, and consequently V is decreasing 
quickly. However, soon the contributions by lighter fields under the slow-roll evolution become 
dominant and there is no more rapid change in H, so V becomes flatter afterwards. Therefore, 
because of this drop-out effect, typically the spectrum is redder than the spectrum of the single 
field inflation where such an effect cannot happen. 



4.2 Marcenko-Pastur law 

The mass of string axion depends on the details of the compactification which involves many 
factors not yet (or formidable to be) computed. Nevertheless, in spite of the unclear microscopic 
physics, the mass distribution is known to follow the Marcenko-Pastur law [8]: the probability 
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Figure 3: Plots of the power spectra V calculated by Eq. (BSj) corresponding to the calculations 
in Fig. [21 We follow the mass distribution given by Eqs. fl30|) and fl32|) in the left and right 
panel, respectively. 

density function for iv? is given by 



(l + v^)^-m2 ni^-m?[l-^y 



m? 



pK) = a-2 2 ' (32) 

where m is the average value of the mass, i.e. (m^) = m^, and /5 is a model dependent parameter 
whose favourable value is expected to be 1/2. Note that for this case masses are more densely 
packed than the exponential distribution we have used in the previous section: we have 

mi = 1.707110m, 

mjv = 0.292893m (33) 

for /3 = 1/2. 

The evolution the of fields according to the Mar cenko- Past ur law given by Eq. ( l32l) is shown 
in the right panel of Fig. [2l Since masses are densely packed, the fields are evolving faster than 
those under the distribution Eq. (1301) we have used in the previous section, where the masses 
spread over two orders of magnitude. The qualitative evolution of the fields is, however, not too 
different: heavy fields exit the inflationary regime at early times, and light fields drive inflation 
at late stages. 

The relative dense packing of the masses is easily read from the power spectrum V and 
the spectral index rig, which are shown in the right panels of Figs. [3] and HI respectively. In 
Fig. [31 we can see that the amplitudes of the spectra are relatively large and fiat compared with 
the corresponding curves in the left panel where the exponential distribution of the previous 
section, Eq. f[5U]) . is applied. Larger amplitude is easily anticipated from the fact that masses 
are densely spaced: the most massive field is 100 times heavier than the lightest one in the 
left panel, while now it is less than 6 times. Hence the energy density is large with the same 
initial conditions thus correspondingly V is also large. Flatter spectrum could be understood 
in the same way: since the masses are very close to each other, the drop-out effect is weak and 
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Figure 4: The spectral indices ris of the spectra shown in Fig. [31 The shape of the spectrum is 
identical irrespective of the overall mass scale. In the left panel, the mass distribution follows 
Eq. (!30|) . with the lightest mass being 1/100 of the heaviest one: for the case of iV = 250 where 
we can obtain A/" ~ 60, rig is decreasing quickly at the early stages due to the "drop-out" effect 
of the massive fields. For the Marcenko-Pastur law given by Eq. ( l32l) where the masses are 
much more densely spaced, shown in the right panel, this effect is very weak and Ug looks more 
or less the same at 60 e-folds before the end of inflation. Still we can see that the index for 
= 250 case is redder than ris for = 500 and = 1000 cases. 

hence the spectra with different numbers of fields look very similar. This is clearly seen by the 
spectral index shown in Fig. HI Unlike the case of the exponential distribution shown in the left 
panel, ng's look very similar now. Still, ng for the case of 250 fields is slightly redder than the 
other cases where more fields support inflation so the effect of the slow-roll is more powerful. 
We have summarised the results of the current and the previous sections in Table |2l 

4.3 Comparison with single field case 

As we have seen in the previous sections, the number of the axion fields and the corresponding 
mass distributions have important effects on the power spectrum and the spectral index. And 
this affects the physical quantities at the end of the inflationary epoch we have calculated in 
Section [3l To get a clearer idea, in this section we compare two mass distributions we have 
discussed in the previous sections for multi-field inflation with the case of single field inflation. 

We have performed numerical calculations on the basis of the single field inflation of Sec- 
tion[2], where we set (pi = 20mpi and m = 10~^mpi. We have tried total 400 fields for multi-field 
cases, and for comparison we have set the average mass the same: for the exponential distribu- 
tion, Eq. ([30]), 

mi = 10m4oo = 2.15313 x 10~^mpi , (34) 

i.e. the lightest mass is 1/10 of the heaviest one. We simply set rh = 10~^mpi for the Marcenko- 
Pastur law, Eq. ( !32l) . then we have 

mi = 1.69594 x 10~^mpi , 
771400 = 0.299126 X lO^^mpi . (35) 
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We have set all the initial amplitudes the same, 



c{0) 



l.OOOOOmpi, 
0.999968mpi 



(36) 



for the exponential distribution and the Marcenko-Pastur distribution respectively, making the 
initial Hubble parameters ifini's for all the three cases identical. In Fig. [5] and Tabled! we show 
the results of the calculations. 
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Figure 5: Plots of the power spectra V (left panel) and the spectral indices (right panel). 
Here we set A/" = to be the beginning of the inflationary phase. At A/" = 0, the amplitudes of 
the power spectra for all the three cases are the same, but the spectral indices are not: this is 
because for the exponential distribution the masses are sparsely placed in heavy region, so the 
drop-out effect is larger at early times. Note that since the masses are densely packed in light 
region for the exponential distribution, more fields support inflation at late stages than for the 
Marcenko-Pastur distribution and hence becomes larger before the end of inflation. 



In the left panel of Fig. O we plot the power spectra V . Since -ffjni's are the same, the 
amplitude of P's are identical initially. However, as inflation proceeds they drop at different 
rates: V for the single field case is most slowly evolving, while for multi-field cases we obtain 
steeper spectra. Here, how quickly the spectrum is decreasing depends on the density of the 
masses: if the mass density is high, i.e. masses of the fields are densely packed, the drop-out 
effect is alleviated and the spectrum is slowly changing compared with the case of sparse mass 
density. This is clearly seen in the right panel of Fig. [5] which shows the spectral indices n^. 
Note that at later stages, around A/" ~ 80, ng under the Marcenko-Pastur law becomes steeper 
than the one under the exponential distribution: following Eq. (|30l) . the masses of the fields 
are densely placed in the lighter region. Thus at later stages where light fields support inflation 
the energy density under the exponential distribution is changing more mildly, making the 
spectrum bluer. 
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5 Conclusions 



In this paper we have examined the dynamics of multi-field infiation at its final stage and the 
density perturbations. Within the slow-roll regime, the fields follow the evolution equations 
Eq. flTTl) according to their masses and we can write all the fields and their derivatives in terms 
of the lightest field as Eq. (fT9|) . Then, the Hubble parameter at the end of infiation, Hend, 
becomes independent of many factors, e.g. initial values and the number of the fields and the 
detail of the mass distribution, and is approximated to the lightest mass tun only, given by 
Eq. (I2T1) . Also the amplitude of the lightest field at the end of infiation, 0^^^, depends on the 
number of fields which have not yet relaxed to the minima of the potentials as Eq. (l26l) . The 
power spectrum V and the spectral index Ug given by Eqs. (l28ll and ( |29l) respectively, however, 
cannot be calculated in the same way because many fields other than the lightest one contribute 
at earlier stages of infiation. 

The detail of the mass distribution greatly affects V and Us'. we have investigated the 
exponential distribution, Eq. (l30i) . and the Marcenko-Pastur law, Eq. ( 132|) . It seems that three 
factors of the mass distribution, namely, the number of the fields A^, the overall mass scale, 
and the density of masses in the range over which they are spread. The amplitude of the power 
spectrum V depends on the overall mass scale, and the spectral index Us is independent of the 
mass scale but is affected by the two other factors. The "drop-out" effect is crucial here: the 
most massive axions fall to their minima at the early stages of inflation and their energies are 
dissipated by the large Hubble parameter, completely dropping out of the inflationary regime 
and contribute no more. This leads to the relatively steep decrease in H and hence makes rig 
redder than the single fleld case. If the masses are densely packed and the number of the flelds 
are large, the drop-out effect is alleviated because soon the slowly rolling flelds dominate the 
energy density. 

Note that, however, detail of the initial conditions and the masses of the axion flelds depend 
on the structure of the microphysics of string compactiflcation. Although it is a formidable 
task to calculate any speciflc and concrete realisation from the flrst principle, this would be an 
important subject to implement a consistent and successful scenario of inflation in the context 
of string theory. 
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Table 2: Summary of the numerical calculations we have performed in Sections 14.11 and 14. 2[ For analytic estimates, we have 
used the results derived in Section [3l Note that the mass given here denotes the heaviest mass, mi, and the average mass fh for 
the exponential distribution and the Marcenko-Pastur law, respectively. Also, 0end and H^nd are given in term of mpi, and V 
and Hs are evaluated at 60 e-folds before the end of inflation. 



